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■ j^j • Abstract 

We study the differentiability of Bessel flow p : x — > pf , where (/?f)t>o is BES x (5) 
£\\ ■ process of dimension 5 > 1 starting from x. For 5 2 we prove the existence of bicontin- 

r^i ■ uous derivatives in P-a.s. sense at x ^ and we study the asymptotic behaviour of the 

1) | derivatives at x = 0. For 1 < <5 < 2 we prove the existence of a modification of Bessel flow 

having derivatives in probability sense at x ^ 0. We study the asymptotic behaviour of 
^ ■ the derivatives at t = tq(x) where tq(x) is the first zero of (p x )t^o- 



1 Introduction 

The regularity of flows of diffusion processes is an important problem related to the stability 
of solutions of SDEs with respect to the initial value. This problem is well-studied when the 
coefficients of the diffusion equation are regular( cf.KuNlTA [7], Protter[9]). Some results 
for the non-Lipschitz case is given in Ren, Zhang 
^ ■ As is well-known, Bessel squared process of dimension 5 > 0, denoted by BESQ 2 ^ (5) , starting 

from x 2 , is the unique strong solution of the following stochastic differential equation: for all 
x ^ and t > 

g ; X* = x 2 + 2 [ ^X*df3 s + St, (1) 

where (3 = (flt)t^o is standard Brownian motion. For x > 7 > with fixed 7, X x = (X x ) t ^ 
is diffusion process with locally Lipschitz coefficients on (0, +00). Moreover, the derivatives of 
diffusion coefficients with respect to initial value are also locally Lipschitz on the same set. It 
gives, using the comparison theorem, that the flow of BESQ^ (S) processes with x > 7, is a 
diffeomorphisme up to explosion time for derivatives, which is 

r ( 7 )=inf{O0:A7 = 0} 

where inf{0} = 00. We remark that in the case S ^ 2 we have P(r (7) = 00) = 1 , and in 
the case 1 < S < 2 we get P(r (7) < 00) = 1. In general, we cannot expect to establish some 

2 

regularity properties after explosion time for the derivatives. But BESQ 31 (S) process is a very 
special case in which this study may be possible. It should be noticed that being particular, 
BESQ processes appear relatively often: it is so for radial part squared of Brownian motion; 
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the laws of some local times for Brownian motion are related to BESQ process, and the same 
is true for some processes related with running maximum of Brownian motion (see Borodin, 
Salminen [2]). Another important example is given by the trace of Wishart process which is 
also BESQ process (see BRUpQ). In this context we should also mention Dunkl process which 
radial part squared is also BESQ process ( see L. Gallardo, M. Yor[I],[5]). 

The same comments can be made for BES :r (5) with 5 > 0. This process is related to 
BESQ* 2 (5) in the following way: for t ^ 




It is well known (see, for instance pU], chapter XI) that for 5 > 1, the BES X (<5) process is the 
solution of the following differential equation: for all x ^ and t > 

fi = x + /3 t + £zH f±-da (2) 

A JO Ps 

where (3 = (Pt)t^o is standard Brownian motion. For 5 = 1 the BES X (5) process satisfies: for 
all x ^ and t > 

p* = x + {3 t + L° t (x) (3) 

where L^(x) is local time of BES X (<5) process at zero. For < 5 < 1 BES :r (5) process verify: 
for all x ^ and t > 

f% = x + Pt+ ^^v.p. f -Js (4) 

/ JO Ps 

where the integral in the right-hand side is understanding is in v. p. sense. 

.2 

The structure of BES X ((5) process is simpler then the one of BESQ X (8) process in a sense 
that the equations (j2J), (JBD, (BJ do not contain a stochastic integral. This is the reason why we 
focuss our study on the flow of BES*^) processes. We remark that some indications related to 
the regularity property of Bessel flow with 5^2 and x > can be found in Hirsch, SongJH]. 

The aim of this paper is to study the regularity property of the flow of BES X (<5) processes 
with 5 > 1. We will distinguish the cases of 5 ^ 2 and 1 < 5 < 2, and inside of them also the 
cases x > and x ^ 0. 

Theorem 1.1. For 5 ^ 2, the flow of BES !C {8) processes has (P-a.s.) derivatives of all orders 
with respect to x for x > which are bicontinuous in (x, t) on the set ]0, +oo[x [0, +oo[. 



To prove Theorem 11.11 we reduce first the problem to the case of x > 7 > 0, then we do 
localisation and we use the classical results. The case x ^ is very different from the case 
x > from point of view of properties and, then also from technical point of view. In the case 
x ^ the mentionned above procedure does not work and we have to use some identity in law 
and some fine asymptotics of Spitzer type to conclude (cf. Spitzer[13J, Messulam P., Yor 
M.[8J). 

Theorem 1.2. For 5 = 2, the flow of BES C {5) processes has (P-a.s.) derivatives of all or- 
ders n at x = (and a fortiori for x > 0). These derivatives are bicontinuous in (x,t) on 
the set [0, +oo[x]0, +00 [. Moreover, for the derivatives of Bessel flow {pt)t^o,x^o the following 
asymptotic relations hold. 
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a) For n ^ 1 uniformly on compacts in t of (0, +00) and P-a.s. 

( d n cf \ 
lim In ( — ) I In x = +00. 

x^O+ \ V 3x n ' ' ) 

b) For n ^ 1 the convergence in law sense holds: 

where T\(f$) is the first passage time of the level 1 for standard Brownian motion, 

c) Uniformly on compacts in t of (0, +00) and P-a.s. 

lim %£M) = I/,,-, 

x^o+ \ dx n ox J 

where C/ n _i = U\(Ui — 1) ■ ■ • (U\ — n + 2) and U\ is a random variable given by [Efy with 
v = \. 

d) For T > e > and < 7 < l/(n — 1) with the same U n _i as in c) 

lim e( sup x n - 1 |^/M|V = E(\U n ^) 

For 5 > 2, i/ie /low of BES jC {5) processes has (P-a.s.) at x = i/ie derivatives only up to 

the order n < n(S) where n(5) = 2 + ( and of all orders for x > 0). These derivatives 

5 — 2 

are bicontinuous in (x,t) on the set [0, +oo[x]0, +00 [. Moreover for n < n(5) we have: 

a') uniformly on compacts in t of (0, +00) and P-a.s. 

( d n p x \ 
lim ( In ( — — - ) / In x ) = n(5) — n, 
x^0+ \ V dx n ' ' J K ' 

and also the property c) and the property d) with < 7 < u/(n — 1) and v = 25 — 3 in [Efy . 

Remark 1. For the regularity at x = we have the following picture. If 5 ^ 3 then the flow has 

only two derivatives in P-a.s. sense and no derivatives of order n > 2 even in probability sense. 

1 1 

If m e N* and 2 H ^ 5 < 2 H , then the flow of BES^Ytf) processes has exactly 2 + m 

m + 1 m 

derivatives in P-a.s. sense. We remark that the regularity of the flow is increasing as 5 j 2, and 
for 6 = 2 the flow is C°°. The asymptotic relations a), b) give us logarithmic asymptotics for 
n-th derivative of pf . The asymptotic relations c), d) characterize the behaviour of the ratio of 
the n-th and the first derivatives in P-a.s. and L 1 - sense. 

If 1 < 5 < 2, then BES X (5) process touches with probability 1 and the results will be 
different from the previous case. To present the results let us denote as before for x > 

70(2) = inf{s ^ : p x s = 0} (5) 

with inf{0} = +00. 



3 



Theorem 1.3. In the case 1 < 5 < 2 and x > there exists a modification p of the Bessel flow 
in the space D(R + '*, C(IR + '*) with the following properties: 

a) p is bicontinuous P-a.s. and has bicontinuous derivatives of all orders on the set 
]0, +oo[x [0, t (x)[. These derivatives coincide with the ones of (p x ATQ r x \) x >o,t>o- 

b) For each (x,t) with x > 0,t > p has derivatives in probability sense only up to the 
order n < n(5) with n(5) = - — - ; which are bicontinuous in probability. Moreover, for 
n < n(5) we have: 

(Qn y 
ln( \J* {x) )/ln(A I = n(5)-n, 



c) For n < n(5) we have 



n— 1 P t q {x) , d P V T {x) \ _ U n -i 



P lim ( o y r _1 T0(X> ! , 

where U n -\ is the same as in Theorem \l.S\ and JJ\ is a random variable given by [Efy with 
is = 5 -25. 

d) For n < n(5) and < 7 < (5 — 25) /(n — 1) with the same U n -\ and U\ as in c) 



to(x) I I Qyn I Qy 1 / x ~i(n-l) 



In the case 1 < 5 < 2 and x ^ 0, t > the mentionned above modification of Bessel flow has 
the same regularity as for x > 0, t > 0. 



Remark 2. For the regularity in probability sense for x > and t > we have the following 
picture. For 1 < 5 ^ 3/2 the flow has only one derivative in probability sense. For m 6 N* 

and 2 < 5 ^ 2 the considered modification has exactly m + 1 derivatives in 

m + 1 m + 2 

probability sense. We remark that the regularity in probability sense is increasing to infinity 

as 5 | 2. The interpetation of the asymptotic relations is the same as in Remark 1. 



2 Regularity of Bessel flow for 5^2 

In the case x > we begin with some rather general Lemmas. 

Lemma 2.1. Let p x be the strong unique (P-a.s.) solution of the equation (Hjj with initial value 
x. If for each 7 > the flow of p x with x > 7 is bicontinuous and has (P-a.s.) bicontinuous in 
(x, t) derivatives of all orders with respect to x on the set ]j, +00 [x [0, +00 [ , then there exists 
an extension of the flow on the set ]0, +00 [x [0, +00 [ having the same properties. 
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Proof A simple patching with respect to 7 proves the result. □ 

Let now 7 > be fixed and x > 7 > 0. To localise the coefficients of the equation we 
take a bicontinuous version of p x (see [10J, p. 362 ). For < e < 7 we put 

r e = inf{t ^ : p] ^ e}, (6) 

and 

r = inf{t ^ : pi = 0}, (7) 

with inf{0} = 00. To simplify the notations and since 7 is fixed, we do not write that r e and r 
depend on 7. 

Lemma 2.2. Suppose for each e > t/iere exists a bicontinuous version of the flow of the 
process (pf Ar )t>o having (P-a.s.) bicontinuous in (x,t) derivatives with respect to x on the set 
]j, +00 [x [0, r e [ then there exists an extension of the flow having the same properties on the set 
] 7 ,+oo[x[0,r[. 



Proof A simple patching with respect to e gives the result. □ 

Proof of Theorem W.W Using the Lemmas 12.11 12.21 and the fact that Pir = 00) = 1, we reduce 
our study to the process (p? Ar Jt^o with x +00 [ and 7 > 0, where r t is defined by (JSJ). By 
comparison theorem (P-a.s.) for all t ^ and x > 7 

For x, y e]7, +00 [ we denote 

/7?/,X PtAT e PtAT t / q\ 

y — x 

From (j2J) we obtain the following linear equation: 

Wr. = 1 " ^ /" T£ ^JS (9) 

and, hence, the solution 

^expj-^f^}. (10) 

To take the limit as y — > x we use bicontinuity of the flow of p x , the fact that on the interval 
[0, t Ar e [ we have a minoration: p y s ^ e, p^' ^ e. By Lebesgue dominating convergence theorem, 
the first derivative of the flow is given by: 

(5-1) f tAT * _L 



*£r.=«P ^/ 7— (11) 



The bicontinuity of the first derivative follows in the same way using bicontinuity of the flow 
p x and the above minoration. 
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Now, using the expression of the first derivative and the arguments mentionned for biconti- 
nuity of the first derivative, we prove, in recurrence way, the existence and bicontinuity of the 
n-th derivative. □ 



To study the existence and bicontinuity of the derivatives of Bessel flow at x = 0, we need 
a renforced scaling property and some asymptotic results. 

Lemma 2.3. // we consider bicontinuous versions of Bessel processes, then for all c > the 
renforced scaling property holds: 



Proof This is clair result for a finite number of due to simple scaling 

property of Bessel process and the uniqueness of the solution of (j2j). Then, the result follows 
by continuity. □ 

Remark 3. As a corollary of this Lemma, the law of a measurable functional of y-{p#t)t^o,x>o 
is the same as the law of the rescaled functional of \^{pt^ c )t^o,x>o) ■ In particular, 



and we note using Cauchy sequence characterisation of P-a.s. convergence, that P-a.s. conver- 
gence for the rescaled processes is equivalent to the same for the original processes. 

Lemma 2.4. If 5 = 2 and x ^ 0, then as t — > +oo, 

" 1 



(lnt) 2 A (Pi) 

where Ti(/3) is the first passage time of level 1 for a standard brownian motion. 
If 5 > 2 and x ^ 0, then as t — ► +oo, 



: / -, — zi:ds 



a.s. 



If 5 = 2 and x ^ 0, then as t — > +oo, 



i r 1 



^f w* ds "~ + °°- 



Proof The first two asymptotics are well-known. For instance, we can find the proof of the first 
one in [S] and, the second and the third ones can be found in [3J. □ 

Lemma 2.5. If 5 = 2 and t > 0, then as x —> 0+, 

1 /"* 1 c 
(lnx) 2 J [pi) 2 
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where Ti(j3) is the first passage time of level 1 for a standard brownian motion. 
If 5 > 2 and t > 0, then as x — > 0+ ; 

1 /"* 1 , a.s. 2 

-as 



lnxj (pf) 2 (2-5Y 
If 5 = 2 and t > 0, then as x — > 0+, 



i r* i 



, a.s. 
rdS — > — OO. 



In Wo (P?) 2 

T/ie mentionned a.s. convergences are uniform in t on compacts of (0, +oo). 

Proof Take x > and consider the integral on [0, £]. Make a change of variables s = x 2 u and 
use scaling property of Lemma 12.31 Lemma 12.41 with x = 1 and the Remark [31 □ 

Lemma 2.6. Lei a ^ and /3 ^ 0. TTie integral 

1 sj -ds (12) 
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zs convergent P-a.s. iff a(5 — 1) + (/3 — 2) (5 — 2) > 0. In particular, it is convergent when 
a + (3 > 2. The mentionned convergence is uniform in t on compacts of (0, +oo). 

Proof We notice that the integrand in (IT21) is positive and, hence, the integral is convergent 
P-a.s. to a finite limit or to +oo. We remark that for s ^ the first derivative of Bessel flow 



By Ito formula we have: 



r i 

So, performing a time change with A s = rfw in considered integral we obtain that: 

io (Pu) 

j-^y ds = J exp \ afj u - — I c/w 

where a = 2 — /3 and 6 = a(5 — 1) + (/3 — 2) (5 — 2) and (3 standard Brownian motion. It remains 
to note that the last integral is convergent iff b > 0. □ 

To prove the existence of the derivatives of Bessel flow at x = we need an explicite 
expression for the derivatives of Bessel flow at x > 0. For this we introduce h x = (hf)t^o with 

h^-^f-^ds. (13) 
7 



Let now X n be a set of multi-indices: 

n 

l n = {I = (zi, i 2 , ■ ■ ■ i n ) -k ^ 0, i 2 ^ 0, • ■ ■ i n ^ 0, ^ ri r = n}. 

r=l 

For g G C n (M) and / G X„, J = i 2 , • ■ ■ i n ), we introduce differential monomials 
as well as differential polynomials 

P n (g) = c iQ^9) (15) 

where cj are real constants. 

Then using the existence of the first derivative and recurrence arguments we obtain that for 
all n ^ 1 

^ = ^^n- 1 (^J (16) 

with Pq(-) = 1 and r e defined by (J6j). By the same reasoning we can prove that for all 1 ^ k < n 

— = |>X x^ ds - (17) 

where Qi(-) are the differential monomials of the type (HM and for / = (ii, z 2 , • • • ik), jk = 
Ylr=i Z V- Taking the limit as e — > and using the fact that P(t = oo) = 1 we obtain the 
needed formulas. These formulas coincide with the ones obtained by replacing of t A r e by t in 
(USD anddHD. 

Now for n ^ 1 we introduce the integrals 

i>qo /\rx\n 

b n = X n / / \> ds (18) 
Jo (p X s) n+2 

which are convergent according to Lemma 12.61 Moreover, using change of variables s = s'x 2 
we establish that 

r-oo (vl\n 



and that the law of b n does'not depend on x. Let B = 1 and let us denote by B n the quantity: 

£n = X>f (19) 
k=l 

We notice that the law of B n does'not depend on x since it is so for the b^. 
Lemma 2.7. For each n ^ 1 there exists a real positive constant c = c(n) such that 

i dx n &r n 1 

v4s a consequence for each n ^ 1 i/iere exists a constant c = c(n) such that 

\x n P n (h x t )\^cB%. 
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Proof The proof is going by induction using previous formulas for the derivatives. For n = 1 
we have 

*J>»= (4 - i)i jf(8*' fa (2o) 

and, hence, 

1^(^)1 - 1) &i- 
So, we see that the claim is true with c = max(l, 5 — 1). Suppose that for 1 ^ m ^ n 

am n x am !, x 
I m— 1 _ l J t I ^ ^rx nm-1 I „m _ t_ I ^ 

We show that the needed relations hold for m — n + 1. Below we will denote by c a generic 
constant. From the formula (1161) with replacing n by n + 1 and (1151) it follows that for the first 
estimation it is sufficient to majorate each Qj(h x ) with / G X n where 

and J = (ii, «2> • " > *n) with J^™ =1 ri r = n. Since _B„ is increasing sequence we obtain from our 
suppositions that 

\x n Qj(/if) | ^ cB^f 3 ■ ■ ■ 5™" cBl 

and it gives the first and third estimations of Lemma. 
We remark that 

dx n+l - ^ (px)i„ +1 +2 

with / = (ii,i 2 , ■ ■ , i n +i), J2r=i ri r = n + l and j n+1 = Ylr=l V 
Then we take in account the formula (fTUj) to obtain 

QM) = (Y s x y n+1 (Pi(K)) 12 ■ ■ ■ (Pn(K)y^. 

Since Ylr=2( r ~l)v=^ + l — jn+i, we have from previous estimations that 

I x n+i-j n+1 Q^ p x> ) 1 ^ c (Y s x y^B[ 2 Bf 3 ■ ■ ■ B^+\ 

Using this estimation and doing the change of variables in the integrals we obtain that 

d n+1 K 



x n+1 1 



Q x n+l 



< r h- B i2 R 2 * 3 • ■ • R ni "+i < rR n+l 



since 1 ^ j n+ \ ^ n + 1 and 6 r ^ for each r. So, we have the second estimation and it proves 
Lemma. □ 



Lemma 2.8. Let t > and x > 0. For all fixed n ^ 1 £/ie sequences of random variables 
d n h x 

( x n — L\ an d (x n P n (h x )) are convei 
ox n 

uniform in t on compacts of (0, +00 



d n h x 

x n — — ) and (x n P n (h x )) are convergent P-a.s. as x — > 0+. T7ie mentionned convergences are 
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d n h x 

Proof We prove in recurrent way that for all n ^ 1 the sequence of random variables {x 11 — — -) 

is convergent P — > 0+ uniformly in t on compacts of (0, +oo). We remark that 
by formulas (1141) . (1151) this convergence gives immediately the same type of convergence for 
{x n P n {h x t )). 

For n = 1 we have 

,| W) = (, - l) x £ fi (,- 1) jf (21) 

By Lemma [2.61 with a = 1 and /3 = 3 we prove that the integral is convergent P-a.s., uniformly 
in t on compacts of (0, +oo). 

Suppose that P-a.s. convergence, uniform in t on compacts of (0, +oo) is valid for {x k — ^ ) 

with 1 ^ k ^ n. Then by f|T5*|) we obtain that the (x h Pk{h x )) are convergent P-a.s., uniformly 
in t on compacts of (0, +oo), as x — > 0+. By formula ( JT7I) we have: 

^ n+1 ^ v t _9ML d 

Q x n + 1 L, 'J ( p? )2+i„ +1 ■ 

Jg_£ n +i 

where / = (zi, z 2 , • • • , E"=i r ^ — ^ + 1 an d j n +i = S™=i V ■ We show that each term 

in the previous sum multiplying by x n+1 is convergent P-a.s. uniformly in t on compacts of 
(0, oo). For this we remark that the term corresponding to I — (ii, • • ■ in+i) i n the previous 
sum times x n+l is equal to 

x n+i t (n z ) j " +i (A(^)) t2 (p 2 (^)) t3 ■ ■ ■ (p n m) in+i dc 

Jo (p x ) 2 +^ 
and that it is equal in law to 

(y s i)^ +1 (p 1 (/ i ^ 2 )) i2 (p 2 (^)) J3 • • • (Pn(K /x2 )y^ 



(pi: 



ds 



where the Pk(h\/ x -i) are equal in law to the x k Pk{h x s ). 

We notice that for 1 ^ k ^ n the (x k Pk{h x )) are uniformly bounded by B k according to 
Lemma |2~T1 and that the law of B\ does'not depend on x. Moreover, since j n +i ^ 1> the integral 

ds 



x\2+j„ +1 



(p 

is converging. So, changing space, we have P-a.s. convergence by Lebesgue dominated conver- 
gence theorem. Then, the final result follows by Lemma 12.31 

The same can be done simultaneously for the expression of (x n+1 - ) and this proves 
P-a.s. convergence of this variable, uniform in t on compacts of (0, +oo). □ 



Lemma 2.9. For n ^ 1 and t > let U n = lim x n P n (h x ). Then U n ^ (P-a.s.). 

x— +0+ 
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Proof Writing the expression for (n+2)-th derivative of p\ from (|T6|) and derivating the same 
expression for (n+l)-th derivative of p\ we get that: for y > and t > 

Pn + l{hl) = ^- y (h y t )P n (h\) + ^(P n (h y t )). 

We notice that for n ^ 1 

lim ^ +1 P n+1 (/^) = f/„+i, lim y"P n (/i?) = C/„, lim = U x 

y-*0+ y^o+ y^o+ dy 

and we prove that 

U n+1 = {U 1 -n)U n . (22) 

For this we take x > and a e]0, 1[ and we integrate the previous equality on the interval 
[ax, x\. 

P n+1 (h y t )dy = [ X ^{h\)P n {h y t )dy + P n (h* t ) - P n (h°*). 

J ax 

Then we estimate each integral, we multiply the result by x n , we take the limit as x — > 0+ and 
we obtain (1221 . 
Finally, we have 

U n+1 = U 1 (U 1 -l)(U 1 -2)---(U 1 -n) (23) 

with 

We show that the random variable JJ\ has a density. For this we make a change of variable 
as in Lemma [2761 with A s = ———-ds to prove that 

JO (Pu) 2 

Ui = (5-1) f exp ((3 U - ( 26 ~ 3 K ) du 



o 



2 

where (3 is Brownian motion. Using Dufresne equality (see for instance [15], p. 95) we have 



U, & ^ (24) 

where Z u follows gamma law T(u, 1) of index v = (25 — 3). 

Since Z7i has a density, P(U\ G N) = and for each n ^ 1, U n ^ (P-a.s.). □ 

Proof of Theorems \1.S\ We have from (T2]) 



Ps 

For x > and t > we put 

ryx _ Pt ~ Pt 
L t — 5 

and we remark that Z% satisfies a linear stochastic equation with the solution given by: 



Z JO PsPs 
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The fact that for all s ^ 0, p x j p° s (P-a.s.) as x j and the property: P-a.s. for alH > 

,t 1 



-ds = +oo, (25) 



/o (p°) 2 

together with Lebesgue monotone convergence theorem give that P-a.s. for all t > 

d -§- \ x=0 = lim 2% = 0. 

OX x^0+ 

In the same manner we establish that P-a.s. for all to > 

lim M = l im exp{-^l /' = 

^ 2 Jo (P x s ) 2 

and this proves the continuity of the first derivative. 

To prove the existence of the n-th derivative at x = equal to zero we show that there 
exists n(5) such that for 2 ^ n < n(5) (P-a.s.) 

lim - V-T- = °- ( 26 ) 

To find n(5) we write that 

an-1 „x 

where h x is defined by f iTSj) and P n is given by (|15[) . From Lemma [2.81 we have that uniformly 
in t on compact sets of (0, +oo) 

lim x n - 2 P n _ 2 (h x ) = U n _ 2 
where £/ n -2 is different from zero with probability 1. This means that ( 1261) is equivalent to 

-(5-i) r i 



lim exp 

2^0+ 



o (p x s y 



-ds + (n — 1) \nx 



(27) 



If 5 = 2 then applying Lemma 12.51 we see that the last relation holds for all n ^ 2 and we can 

put n(S) = +oo. If 5 > 2 then it is easy to see that for n(8) = 1 + - — - and n < n(S) the 

o — 2 

relation ( 1271) holds and for n > n(S) it fails. If 5 > 2 and n = n(S) then by Ito formula 

ln(pH = lnx + / — d/3 s + / T1 ^ds. 

Jo P x s 2 J (pf) 2 

We apply a central limit theorem (see |12j , p. 472 ) for the martingale for M = (M t )t^o with 

1 f t/x2 1 

obtained from original one by time change, to prove via Skorohod representation theorem that 
the quantity appearing as the power in exponential in ( 12T1) . namely 
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behaves as c£| lnxj 1 / 2 as x — » where £ is standard jV(0, 1) random variable and c is some 
positive constant. Hence, (j2*7|) fails on the set {£ > 0} of probability 1/2, as well as (12U|) . 

For the bicontinuity of the n-th derivative at x = for 2 $C n < n(5) we prove that P-a.s. 
for all t > 

hm ^4 = 0. (28) 

t—*to 

The proof of (j28j) is going in the same way as f)26p using the fact that the convergences in 
Lemmas 12.81 and 12.51 are uniform in t on compact sets of (0, +oo). 

The asymptotic relations a), a') b), c) follows from (fl6l) and Lemmas 12.51 and 12.91 To prove 
d) we remark that according to Lemma 12.71 we have 

sup \x n - l P n ^{h x t )\^cBlz\ 

with B n defined by f)19p . It remains only to show that B^Lx is integrable. Since for a, b G K + 
and 7 > 0, (a + b) J ^ c(a 7 + W) with some constant c, B^Ti is integrable if }fj} n ~ l ^ k are 

r* i 

integrables for 1 ^ k ^ n — 1. Making time change with A t = k 2 / tig and using Dufresne 

Jo [plr 

identity (see [2], p. 78) we obtain that 

b * = 13. e MPu 777 ) du 



k 2 Jo " 2k kZ v{k) 

where (/3 u ) u ^o is standard Brownian motion and Z v ^) is the variable following gamma law 
T(u(k), 1) of index v(h) = So, we have needed integrability if for all 1 ^ k ^ n — 1, the 

variables (2/Z v ^) * are integrables. As well-known this is true, if 7 ^ ' < v{k) and the 

25-3 
n-1 



last condition is satisfied for 7 < ^— p. □ 



3 Regularity of Bessel flow for 1 < 8 < 2 

We start with some lemmas needed to prove the theorem 11.31 Let tq(x) be defined by (j5J). 

Lemma 3.1. Let x > be fixed. Then tq(x—) = tq(x) = tq(x+) (P-a.s.). Moreover, there 
exists a cadlag version of (tq(x)) x> q. 

Proof By comparison theorem we have that for x ^ y (P-a.s.) 

To(x) < r (y). 

Then, there exist the limits (P-a.s.): 

lim tq(x) = r (x+), lim r (x) = r (x—). 

y-*x+ y^x- 

We take y < x < z then by comparison theorem again for 7 > 

E(r (yy) < E(r (xy) < E^zf) 
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Since 

ro(x) ± (29) 

where 7„ is random variable of gamma law T(v, 1) with index v — 1 — |, we have for small 
7 > that 

E(r (x)y = c 7 x 27 

with some positive constant c 7 . The mentionned informations implies that for small 7 > 

E(r (x-y) = E(r (xy) = E(r (x+y) 

and, hence, P-a.s. r (x— ) = tq(x) = t (x+). 

To construct a cadlag version of (to(x)) x >o, we take for x e Q + the value of tq(x) and for 
x G R + \ Q + we put: 

To(x-) = lim r (y), r (x) = lim r (y). 

y—>x — J/^a: + 

J,eQ+ J/GQ+ 

This construction gives a cadlag version of (r (x)) :E>0 which preserves the finite-dimensional 
distributions of t (x). □ 

Lemma 3.2. We consider bicontinuous (P-a.s.) modifications of {pt)t^o,x>o- Then the flow 
of processes (Pt/\ TO ( x ))t^o is bicontinuous in probability and has the first derivative in probability 
sense with respect tox. This derivative is bicontinuous in probability on the set]0, +oo[x]0, +00]. 



Proof Let (x, t) be fixed with x > 0, t ^ 0. To prove bicontinuity in probability of (Pt/\ To ( x ))x>o,t^o, 
we write: 

\PsAT (y) ~~ /W (si) I ^ \PsAT (y) ~ PsAT {y) \ + \PsAT {y) ~ PtAT (x)\- 

We introduce the set Ag = {\t(x) — r(y)\ ^ 5} with < 5 < 1. Then, on the set Ag for 
\s — t\ < 7 < 5 we have: 

\P V sAT (y) ~ PtAro(x)\ < SU P \PI ~ P X u\ + S U P \Pu ~ Pv\- 

O^u^t+S \u— v\^S; u,v^t+5 

We write for e > that 

n\p y S AMv) - p^a p({\p y SAMy) - P x t A Mx) \ >€}n As) + p{ai). 

From the previous estimations we obtain that 

P (\p y sAM y) -PtAMx)\ >e)^P( sup \p\-pl\ > c/2)+P( sup \pl-pl\ > e/2)+P(A c s ). 

0<u<t+l \u-v\^S;u,v^t+l 

Using the facts that (p^)«^o,x>o is bicontinuous P-a.s. and that t(x) is continuous in probability, 
we obtain taking lim^o lim y ^ x , the claimed bicontinuity. 

We show that for t G]0, +00] the first derivative of (p^ Ar ( x ))t>ox>o with respect to x is given 

by: 

(S — l) r tAT °( x ) r/ s 
Y t U x) = eM-^ I ^} (30) 
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For this we take e > and we do a localisation with 

T e (x) = inf{s ^ : p x s ^ e}. 
Then we write ( ITU]) . We notice that r (x) = linv_> r e (x) and, then, 

^.^-Pf-V 1 / ^?>- (31) 

z JO PsP s 

where Z x ' y is defined by ([H])- Via comparison theorem it can be shown that 

: tAr (x) dg f tAr (x) ^ 

lim 



v-^Jo PsPs Jo (Ps) 2 

and taking lim,,^ in (I3T1) we have (IHUl) . To obtain the result for t = oo it is sufficient to 
take lim in (jTUl) and continue in above way. 

By time reversal we have: 

£ (Kij-flo^^fx)) = £ (iPt A ~ S )o^L(x)J (32) 

where p x ' s is BES Z (<5) process and L(x) = sup{t ^ : Pi' 4 " 5 = x}. Since 1 < 5 < 2, we have 
that 4 — 5 > 2 and using asymptotics of Lemma [231 we obtain: (P-a.s.) 

teP = +0 ° (33) 

and it gives together with (13"U|) the expression 



^tAT (o;) 



ex P {-^C W ^} if t<r (x), 
if t ^ r (x). 



Now we prove a bicontinuity of the first derivative at each point (x,t) with x > 0,t > 0. 
We consider tree sets D±, D 2 , D$: 

D 1 = {lu: t (x) > t}, D 2 = {lu: t (x) < t}, D 3 = {uo : r (x) = t} (34) 

and we prove a bicontinuity on each of them. For D\ we write that D\ = [J e>0 D\ where 
D\ = {r t (x) > t}. On each set D\ bicontinuity of ^Ws) f° nows from bicontinuity (P-a.s.) of 
(Pt)t>o,x>o- Hence, taking a countable set of e we obtain the same result on D\. On the set D 2 
we have ^Wz) = ^toO) = ^' anc ^' hence it is continuous. 

Take now the set D 3 , then t = Tq(x). Let (s, y) be in the neighbourhood of (r (x),x). 
We show that YsAT (y) * s * n ^ ne neighbourhood of = 0. In fact, on the set {s ^ r (?/)}, 

YgAroiy) = ^rofe) = ^' ^ n ^ ne se ^ ^ s < r °(^)} we remar k that for all x > (P-a.s.) 



y—>x 



I du 
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In fact, by comparaison theorem for small 7 > and e > we have 

.^iw'i"! w = l w (36) 



Taking e — ► and then 7 — > we have from flHoj) and (136]) the relation (|35[) . □ 

To investigate the existence of the derivatives of higher order, we prove the following two 
lemmas. 



Lemma 3.3. Let x > be fixed. Then 

y^ x (y — x ) 



plim ro(y)-r (x) =a 



Proof Let y > x > 0. Since (p^( x \, z > 0) is J r T0 ( x ) - Markov, we have 

T (y) - r (x) = mf{u > : = r (^ o(x) ) (37) 



Using fl29l) we obtain that 

where 7^ is gamma random variable with index u = 1 — 5/2, independent from p y TQ ^ x y But 

(^tqO)) 2 _ (Pr (x) ~ Pt {x)) 2 



(y - x) 2 (y - x) 2 

since P* Q t x } = 0. It was shown in Lemma [3.21 that (P-a.s.) 



(39) 



lim ^ ^ = 0. 



y-+a: y — £ 

Then, ()37|). (1551) and (1591 implies 



2/->a;+ (?/ — x) 2 



The same consideration with x > y > gives again (159]) with the exchanging x and But 

i-y *=<»> F1 2 J„ 

and the relation (P — a.s.) 

r {y) ds 

lim / — -y -> +00 



implies 

and it proves the result. □ 



P lim T ° ( f } ~ T j X) = 

(y — x) 2 
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Lemma 3.4. The flow of the processes (^ TO (x))*>o defined by K3(J\) has bicontinuous derivatives 
in probability sense at x > only up to the order n < n(5) where n(5) = - — -. 

Proof First of all we remark that D\ = [J e>0 D\ and on the sets D\ the flow of {Y t x ATQ ^)t>o with 
x > has the bicontinuous derivatives of all orders. This follows from the fact that on D\ this 
process coincide with (Y^ T j x \)t>o,x>o an d we can use the previous results for classical case. On 
the set D 2 the result is also trivially true. 

Let (x, t) E D 3 be fixed with x > and t = r (x). Since Y^,-^ = we evidently have that 

QnyV 

p lim ^(y) = Q 
y-*x- dy n 

If we show that there exists n(S) > such that for n < n(5) 

P i im * Ar ° fa) = 0, (40) 
y^x+ dy n 

then the mentionned relations and continuity of the derivatives on D\ and D 2 will imply that 
the flow of the processes (X t x AT0 ^ x ))t>o has continuous derivatives of the order n < n(5). We 
recall that for y > x and t = tq(x) 

QnyV 

= y^K^) (41) 

where h v and P n are defined by (jT3"j) . (|ToT) . 

Let u = y/x. Performing time change s = s'x 2 we obtain 

Ki») p »Kw) = WnKw)/^ (42) 

First of all we investigate the behaviour of 



Using time reversal we obtain that 



r (l) i \ / pL u {v) i 

^1/ 7^**\ptto=A=£\ T^-Js\ (43) 



where L u {v) = sup{s ^ : p v s ' A ~ 5 = u}. By time change we obtain that 

L u(v) j j-L u/v (l) ^ 



o (p. ) 2 (p s ' ) 2 

where L u /„(1) is defined as previously with replacing of u by u/v and t> by 1. Since as a — > +oo 

-^j--*>(l) 
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where t (1) is the corresponding time of attending of zero, we have that 

hiL a (l) , 



P lim 

a^+oo 2 In a 



Then using Lemma 12.41 we obtain that 

P lim , . ., / . , , ds = - -. (45) 

Since lim p" m = 0, we obtain using standard arguments from (j45j), (T43l) and (jHJ) that 



1 r^ 1 ) 1 2 

P lim — / ^^ds = -. (46) 

u-+i + \ln(pl (1) )\J (p«)2 2-6 y } 



If we appy a time reversal to P n (h™ o ^) then we obtain 
From Lemma 12.81 and 12.91 we obtain that 



v,4-8 
(v) 



P lim v n PJh V r A 7i) = U, 



where U n is defined by formula (1231) with XJ\ given by ff24|) and = 5 — 26. So, by standard 
arguments we deduce that 

P Km + (pl (1) r P n (hl (1) ) = U n (47) 
Since U n ^ with probability 1, we obtain finally that the relation (|4T?|) is equivalent to 

"(5-i) r o(1) i 



lim exp 



I 2 U ^ dS + nHplilh 



0. (48) 



Let n(5) = (6 - l)/(2 - 6). It is easy to see from g6]) and g?D that if n < n(6) then (I15T) 
holds and that if n > n(6) then ( 1481) fails. For n = n(6) it can be shown using a central limit 
theorem for martingales and Skorohod representation theorem like in the proof of theorems 11.21 
that the limit in (1481) exists only on the set of probability 1/2, and, that (1481) fails, too. 

□ 

Proof of theorem \1.3\ Let p° be continuous version of BES°(5) process starting from zero and 
(to(x)) x> o be a cadlag version of the corresponding process. For all x > and t ^ we put: 

Pt = PtAT (x)h0Mx)i + (Pt ~ Pro(xo) J N(x),+oo[- (49) 

Using strong Markov property we prove that the both processes p and p have the same finite- 
dimensional distributions. 

We show that the trajectories of (I49p are in P(R + '*, C(R + '*)). For this we remark that for 
each x > 0, the process (Pt)t^o is continuous in t. Moreover, for y > x 

\p y t-pt\^2 sup |p^-p*| + sup \p° u -p° v \ 

u,v£U(x,y) u,vGU(x,y) 
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where U(x,y) = {(u, v) : \u — v\ ^ (to(v) ~ T o( x )) ', u,v ^ t + (r (y) — r (x))}. Since for all 
x > , 

lim r (y) = t (x+) = t (x), 

we see, that uniformly on compact sets of t, the right-hand side of the last inequality is tending 
to zero as y — > x+. Taking y < x and using the fact that for all x > 

lim T (y) = to(x-), 

we obtain the existence of left-hand limits uniformly on compact set of t. 

We show that the first derivative of p coinside with the one of (pf Aro ( x ))t^o,a:>o- I n fact, 
consider tree sets Di,D 2 ,D 3 defined by fl34l) . On the set Di the process (pf)t^o coincide with 
(PtAT (x))t>o,x>Oi an d the existence of the first derivative was already discussed in Lemmas 13. 2\ 
13.41 On D 2 the same process coincide with (pt)t^ TO (x),x>o according to comparison theorem and 
the first derivative is equal to zero. On the set D 3 we have t = t (x) and 

y x 

^ -7p — — lf y > X, 

rt Pt_ _ J y — X 



y-x 



We obtain as in Lemma [3.31 that 



if y < x. 



lim 



y x 

Pt (x) - Pt (x) 



y^x+ y — X 

Hence, the first derivative of p coinside with the one of (Pt^ TO ( x ))t^o,x>o an d we obtain the claims 
from Lemma 13.41 and from the relations (jHJ) , ( )42l) , (HHJ) and ( 1471) . 

Now we consider the case x = and t > 0. We put (p°)t>o — (Pt)t^o and we remark that 
for x > 

Pt ~ Pt = PtAro(x) ~ PtAro(x) (50) 

Making time change and using scaling Lemma we have: 

i x o \ £ l o 

~\PtAr (x) ~ PtAro(x)) - P(t/x2)Ar (l) ~ P(t/x 2 )AT (l) 

Since P(r (l) < oo) = 1, we obtain that for t > 

P l im PtlA = o. 

x^0+ X 

Then we verify easily a bicontinuity of the first derivative at x = and t > 0. For x > we 
have 

v x — v x — v 1 

I t — I tAr (x) — I t/x 2 Ar (l)- 

Since P(ro(l) < oo) = 1 and ( |33l we obtain that for each e > 

P(Y t %>AMD > e ) = P ( 7 "o( 1 ) > l l x ^ ^ 

as x — > 0+. 



19 



From the formula (fl~6j) we have for x > and 1 ^ n < n(S) that 

Pt_ _ ^Arojx) _ vx p (hx . 

dx n dx n ~ l I tAT (x) r n-l[.'l'tATo(x)) 

and, hence, 

~~Q^n ~ Y (t/x^Ar (l)^n-l[H {t/x2)ATo{1) )/X . 

We can see from the Lemma [2.61 and the proof of Lemma [2.81 that for 5 > 3/2, Pn-i(h^m) is 
a finite random variable. So, for each e > 

p(-^A > c) = P(To(1) > t/^ ^ o 

as a; — > 0+. It means that for 3/2 < 5 < 2 there exist the bicontinuous derivatives of order 
1 ^ n < n(5) with respect to x in probability sense at x ^ and t > 0. □ 
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